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Abstract 

We study the effect of advection and small diffusion on passive tracers. The ad- 
vecting velocity field is assumed to have mean zero and to possess time-periodic 
stream lines. Using a canonical transform to action-angle variables followed by a 
Lie-transform, we derive an averaged equation describing the effective motion of 
the tracers. An estimate for the time validity of the first-order approximation is 
established. For particular cases of a regularized vortical flow we present explicit 
formulas for the coefficients of the averaged equation both at first and at second 
order. Numerical simulations indicate that the validity of the above first-order esti- 
mate extends to the second order. 
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1 Introduction 



The characterization of the behavior of a passive but diffusing scalar advected 
by a prescribed, smooth velocity field has been the subject of intensive research 
going back at least as far as Batchclor [3]. Above and beyond the obvious prac- 
tical importance in applications ranging from micro-mixers to global climate 
dynamics, 'scalar turbulence' as exhibited by solutions of the linear advection 
diffusion equation also provides an avenue for insight into the structure of the 
Navier-Stokes equations [20] . 
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Despite the linearity of the governing equation, complete characterization of 
scalar solutions, especially the asymptotic decay of such solutions for vanish- 
ing diffusivity, continues to pose considerable difficulties even when restricted 
to the case of planar flows. Here we are concerned with the so-called Batche- 
lor regime, where the spatial scale of the velocity (ly) is assumed to be much 
larger than the diffusive length scale (l^)- For the scalar turbulence problem 
with periodic velocity fluctuations in time and space, rigorous homogeniza- 
tion techniques [4, 15, 17] can be applied to compute effective, renormalized 
diffusivities on large time and space (L 3> Iv) scales assuming that the ini- 
tial distribution of the scalar held satisfles the scale separation Ig <^ ly On 
the other hand, the situation for scalar flelds with variations commensurate 
with both the velocity fleld and the domain size, {Ig ^ Iv L ^ i^) requires 
different techniques. 

The present paper is motivated in part by the phenomenon of persistent pat- 
terns, termed 'strange eigenmodes' [18], that occur under the action of periodic 
stirring. Such patterns, characterized by exponential decay of the scalar vari- 
ance and self-similar evolution of scalar density functions, have been observed 
both numerically and experimentally [22], [5]. Theoretical predictions of the 
decay rates of the scalar and the connection between the observed eigenmodes 
and the phase space of the underlying advection dynamics have been investi- 
gated [1,7,10,19,21], most often in the context of non-linear maps. 

In the strange eigenmode regime, the decay of the scalar contrast can be stud- 
ied via Floquet theory. While well established for ODE's, the existing theory 
for parabolic PDE's requires that the PDE satisfies a restrictive spectral gap 
condition which the advection diffusion equation fails for vanishing diffusivity 
(see [14]). An alternate approach, shown by Chow et al. [6] for one dimensional 
parabolic equations, is to prove the existence of an inertial manifold for the 
system and then apply ordinary Floquet theory to the inertial form. Our goal 
here is less ambitious, but a potential first step in this direction. Following 
Krol [13], we propose a formal averaging procedure for the advection-diffusion 
equation when the velocity field is has zero mean and possesses time-periodic 
stream lines. The approach is perturbative, making explicit use of the dispar- 
ity of time-scales between the advective and diffusive operators. Advection 
fields of the form u = u{^,t) = u{$,)f{t) guarantee that, in the case of van- 
ishing diffusion, the time-dependent system can be solved using action-angle 
variables, and that tracer trajectories will be time-periodic. The explicit solu- 
tion in action-angle coordinates allows the original equation to be written in 
a form suitable for averaging. By applying Lie transform techniques, we de- 
rive an approximate averaged equation. The fact that this resulting equation 
- in contrast to the original problem - has time-independent coefficients, facil- 
itates the theoretical and numerical analysis tremendously. The use of the Lie 
transform also allows relatively straight-forward computation of higher order 
corrections. 
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The form of the paper is as follows. In section 2 we consider the transformation 
which places the advection-diffusion equation in a form suitable for averaging. 
Lie transform techniques are used to average the equation in section 3. A 
proof of the convergence of solutions of the averaged equations to those of 
the original time dependent problem is given in section 4. An application to 
a specific flow field, a periodically modulated, regularized vortex, along with 
numerical comparisons of the solutions are given in sections 5 and 6. 



2 Action-angle variables 



We consider the advection diffusion equation in the following form 

Q + (m ■ V)c- kV^c = 0. (1) 

All functions depend on the spatial variable ^ = {x,yY and the time t. We 
look at (1) as an initial value problem, assuming that 

c(0,O=c(^nO (2) 

is a known function. Incompressibility V-m = implies that the given velocity 
field u is derived from a stream function \E' such that 

w(e,t) = v^vi>(e,t) (3) 

where V"*" = {dy, —dx). 

We assume that the stream function \1/ is of the particular form 

^{U) = ^iOfit) (4) 
where the function / is periodic in time with period T. For consistency in the 



averaging which follows, we also require that (/) = ^ f{t)dt = 



A standard non-dimensionalization of (1) with velocity, length and time scales 
given respectively by (f/, L, T) gives 

Q + -^(M-V)c-eV2c = (5) 

where the dimensionless groups are the Strouhal number, St = UT/L, the 
ratio of the forcing period to the advective time-scale, and e = Tn/L"^, the 
ratio of the forcing period to the diffusive time-scale. Throughout, we assume 
St = C(l) and e -C 1. To clarify notation, we take u = u/St throughout. 

Our aim is to derive an equation of the form 

Ct + Cc = (6) 
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with a time-independent local linear operator £ that can be used to construct 
approximative solutions to the initial value problem given by (2) or (5), in 
the limit of small diffusivity k -C 1. Due to the fact that the leading order 
evolution is given by the periodically varying advection operator, we cannot 
apply averaging techniques directly. Instead, we seek a transformation to the 
Lagrangian frame which results in a new equation where the coefficients of 
both the advective and diffusive operators are zero-mean periodic functions of 
time. The transformed equation is then in a form suitable for averaging. 

For the restricted class of flow fields considered, the the proper transformation 
is simply to action-angle coordinates of the underlying, conservative advection 
equation. Introduce the function F as 



and write the tracer coordinate ^ 
autonomous Hamiltonian system 

dF 



= f fit') dt' (7) 



= ^^^iO ■ (8) 



function of F. We then obtain the 



Since this system integrable, there exists a canonical transformation [2] 

C:{x,y)-^{J,9) (9) 

such that the advection-diffusion equation (1) becomes 

ct - f{t)uj{J)ce - e (r : VV + (5 ■ V) c = (10) 

with a matrix F = r(^^, J) and a vector 6 = 6{6, J), that are solely determined 
by the canonical transformation C. The advantage of the representation (10) 
lies in the fact that the evolution of the unperturbed problem is linear and 
given by 

j = Jo, e = eo-uj{j)F{t). (11) 

Therefore, we can now use these stream lines as coordinates via the transfor- 
mation 

c{t, j,e) = v(t,j,e = e- cu( j)F(t)) (12) 

and with the transformation rules 
cj = -u'Fvq + vj, 

cjj = {uj' fF^Vgg - 2uj'Fv^j - uj"Fv-Q + vjj 
and the rescaling of time as r = et the equation for v takes the form 

Vr = {V -.VW + ~5 ■V)v . (13) 
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Since the only explicit time dependence in the coefficients F and 5 is given 
in terms of the T-periodic function F, the equation (13) is now suitable for 
averaging. 



3 Lie transform averaging 



In order to average (13), we use a technique based on Lie transforms ffist 
developed in the finite-dimensional context [16] and then applied to cases 
involving an infinite number of degrees of freedom [9, 12]. The basic idea of a 
Lie transform is to use a near identity transform of the type 

v = exp(0 ■ Vl)\/ . (14) 

The linear operator, 0- Vl, is chosen to eliminate the explicit time dependence 
of the coefficient of an equation 

Vr = X{v,T) (15) 

in order to obtain an equation with time-independent coefficients of the form 

Vr = Y{V) . (16) 

Since the functionals X and Y depend on v and all its spatial derivatives, the 
operator ■ will be defined in our case as 

(t)-^L = ^(t)nx,my-^. (17) 

n,m (yVnx,my 

where (j)nx,my = d^'^^'^^cj)/ dx'^dy'^ and Vnx,my = d^'^^'^^V/dx'^dy'^ respectively. 
The subscript at Vl distinguishes this operator from the usual V. The gener- 
ating function also depends on V and all its derivatives. The idea is that the 
explicit time dependence will be kept in (p rather than in the equation for V, 
hence will also depend periodically on r. The general transformation rule 
under which (15) transforms to (16) using (14) is [11] 

Y-V+ (£e^'^^) e-<^-^ = e-^-^^ {X ■ VL)e-'^-^^ (18) 

Both terms can be conveniently expanded using the Campbell-Baker-Hausdorff 
formulae 

(^e<^-^-) = (0, + i[0,0,]^ + i[0, [0,0,]z.]l + ...) ■ V, (19) 
e'P-^L^X ■ Vi)e-<^-^- = + [0,X]i + ^[0, [<P,X]l]l + ■ Vl (20) 
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where the Lie commutator is defined through [A, i?]^ = {A ■ Vl)B — {B ■ 
Vl)A and again the subscript distinguishes the Lie commutator from the 
usual commutator. We now expand both Y and in a series in the small 
parameter e as 

Y = Y, + Y^ + ..., = 0, + 02 + ... (21) 

where 0{Yn) = C'(0n) = e"" and difi^erentiation by r lowers the order of 0„ by 
one 

0{d4>n/dT) = e"-i . (22) 

The equation for Y can then be solved order by order. At the leading order, 
we find 

yo + ^ = X (23) 
and averaging this equation yields directly Yq = {X) due to the periodicity of 

01. 

The transformed advection-diffusion equation (13) can be written in the form 

Vr = Lv, L = f : V V + 5 ■ V . (24) 

Averaging this equation immediately yields at the leading order 

Vr = {L)V. (25) 

Here, and in what follows, (...) denotes averaging over one period. The aver- 
aged equation is then 

K= ((f) : VV+(5)-V)y. (26) 

Thus, the time- dependent coefficient are simply replaced by their time aver- 
ages. The leading order of the generating function 0i is found as 



0i = Li\/= ^ L-{L)^V. (27) 
Introducing Fi and as 

f^f-(f). f^S-iS) (28) 

we can write Li explicitly as 

Li = Fi : VV + (5i ■ V . (29) 

Higher order corrections can be calculated in an elegant way using the Campbell- 
Baker-Hausdorff formulae. For the second order term in the expansion of Y, 
for example, we find 
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Y,= ^{[L,V,L,^V]l) + [{L,)V, {L)V]l 
= 1{LL,)-{U){L))V 



(30) 



where the last equahty follows after using the definition of the Lie commuta- 
tor and integration by parts. Collecting first and second order, we obtain as 
averaged equation for V 

K={{L) + e{{LL,)-{L,){L)))v (31) 



4 An averaging theorem for parabolic differential equations 

In the previous section, we applied a technique based on Lie transforms to 
average the equation (13), and we arrived at the equation (25). Here, we state 
and prove rigorously a theorem on averaging of parabolic partial differential 
equations which is due to Krol [13]. We assume that the differential operators 
in (13) are given by 

and 

6 = e{bi{t,x,y),b2{t,x,y)), 

where aij, hi G C^i^B? x [0, oo)), and [aij] is symmetric and uniformly positive 
definite, i.e, there exists ^ > such that for all ^ = (^i, ^2) £ and (x, y, t) G 
R"^ X [0, 00) we have 

2 

a,j{x,y,m,^j>e\^\''. 

In the following, let tq = 0{l/e), and let || ■ ||oo denote the usual supremum 
norm on either or R"^ x [0, tq], depending on the context. 

Theorem 4.1 Letv and V be solutions to the Cauchy problems v{0) = V{0) = 
Vq G C^i^R^) for the equations (13) and (25), respectively. Then 

\\v-VU = 0{t). 

Proof: First note that the existence and the uniqueness of bounded solu- 
tions f, V on C'^{R^ X [0,ro]) is well estabhshed (see [8]). Also, since (25) 
is autonomous, the derivatives Va of V also satisfy an autonomous parabolic 
differential equation with the same second order differential operator, however 
with different smooth and bounded first and zero order coefficients: 

V^r = ((f) : VV + t{br{a){x,y)M{a){x,y)) ■ V + e/(«)(x, y)) V^. 
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The Phragmen-Lindelof principle for parabolic partial differential equations 
implies that 

IlKlloo < IIMalloo + roe||/(a)||oo = 0{1) . 
Let us now define a near-identity transformation 



V{x,y,T) = V{x,y,T) + 



CHs) - (L)) ds 



V{x,y,T). 



Since the integrand is T-periodic with zero average, the equation reads actually 



V{x,y,T) = V{x,y,T) + 



t/T]T 



(Lis) - (L)) ds 



V{x,y,T), 



and it is an easy observation that — y||oo = C(e)- 

On the other hand, one easily verifies that V satisfies the equation 

K = L(r)V + M{r)V, 

where M(r) = Jq{L{s) — {L)){L) — L{t){L{s) — (L)) ds, and the initial-value 
condition V^(0) = vq. Notice that M is a T-periodic fourth order operator with 
smooth bounded coefficients of order C(e^). Consequently, the difference V — v 
satisfies the equation 

{V - v)r = L{r){V -v) + M{t)V, 

and the initial condition {V — v){0) ~ 0. By the Phragmen-Lindelof prin- 
ciple for parabolic partial differential equations, we conclude — t'||oo ^ 
||M(-)\/||ooTo = 0{e). This concludes the proof. 

We remark that similar techniques can be applied to solutions of the second 
order Lie-averaged equations leading to 0{e^) error estimates. 



5 Reguleirized vortical flow field 



To illustrate the theory and to give a comparison to numerical simulations, we 
consider the particular case of a regularized vortical flow field whose stream 
function is given by 

X, y) = In (^a2 + x2 + y2^ f{t) . (32) 

This flow represents perhaps the simplest example for studying the interplay 
between diffusion and nonlinear, time-periodic advection. 
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Since = + is a constant of motion, the unperturbed stream lines are 
given in Cartesian coordinates as 



x{t) = cos {u{r)F(t)) Xo + sin {u{r)F(t)) yo 
y(t) = - sin {u{r)F{t)) xq + cos (u;(r)F(t)) yo 

where 

u;ir) = . (33) 

Obviously, we can take as a canonical transform to action-angle variables 
the usual transformation to polar coordinates (x, y) (r, 6) and in these 
coordinates, the advection-diffusion equation (1) is written for this particular 
flow field as 

ct - uj{r)f{t)c0 = e (^Cr + €„ + ^c^e^ = eAc. (34) 
The transformed equation (13) becomes (using t instead of the rescaled r) 

vt = e (^Av + F (^1^^ + uj"^ + 2u'v-,)j + F\u'fv-,^ (35) 

Using the previous results, we obtain at leading order 

V, = eUv+{F)({^+ A V~, + 2u;%\ + {F'){uJ')X^ (36) 



As shown in the Appendix, the leading order contributions in Cartesian coor- 
dinates produce a time independent advection field with spatially dependent 
rotation and source-like terms. The averaged diffusivity tensor is usually full 
and symmetry breaking. The relative importance of the symmetry breaking 
terms depends upon the explicit form of the time dependence through the 
ratio of (F) and (F^). 

For this flow-field, it is not difficult to compute corrections at second order 
as well. In order to make our notation more efficient, we introduce the two 
operators £i and £2 



+ uj"\dg + 2uj'd-Qdr (37) 




(38) 

and introduce Fi = F and F2 = F"^. The first order averaged equation (36) 
becomes then 

Vt = eiAV+ {F,)CiV + {F2)C2V) 
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Applying now (28), we introduce the functions Gi and G2 that can be found 
exphcitly from Fi and F2 as 

G,(to) = rV.(r) - {F,))dT, J = 1, 2 . (39) 
Jo 

At second order in e the equation is 

Vt = e{AV+ {F,)C,V + {F2)C2V) + ([A, {G,)C,]V + [A, (Ga)^] V 
+((FiGi) - (Fi)(G'i))/:?l^+ ((FsGa) - {F2){G2))ClV 
+((FiG2) - (F2)(G'i))/:i/:2V + mGi) - {F,){G2))C2CiV) (40) 

where [A, B] = AB — BA denotes the usual commutator. 
6 Numerical Simulations 

We can now integrate both the original problem (1) and the averaged equation 
in first order (36) numerically and compare their solutions. In our numerical 
simulations, we use a standard Adams-Bashforth-Moulton Method in Fourier 
space, the Fourier transformations are done using FFTW. We work back 
in Cartesian space and for this purpose, we compare solutions at Poincare- 
sections where F is zero, 9 = 6 and (36) is written in Cartesian coordinates. 
The explicit form of (36) in Cartesian coordinates is given in the appendix. 
Throughout, we consider a single initial condition of the form 

c^''\x,y) = X exp(-6r^) (41) 

where the constants a and b were chosen as a = 1.0 and b = 4.0. 

As shown in (36), the form of the averaging depends explicitly on the nature of 
fit). For the particular case of f(t) = sin(t), we find (F) = 1 and (F^) = 3/2. 
Fig. 1 shows both the initial condition and the evolution after ten periods for 
this choice of f(t). 

In the absence of the time-dependent advection field, the tracer field will obey 
the purely diffusive equation 

ct"^ = eAc^"''^ , c("'^) (0,0 = c(^) (0 (42) 

and the diffusion will simply spread the initial distribution out and preserve its 
symmetries. In the presence of the time-dependent vortical field, however, the 
particles will move back and forth within one period and the interplay of the 
time-dependent trajectories with the diffusion will give rise to a breaking of the 
symmetry of the initial distribution resulting in a "twist" that can be clearly 
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Fig. 1. Evolution of tracers in the time-periodically modulated vortical flow field. 
The figure on the left shows the initial condition, the figure on the right shows the 
result of the tracer motion after 10 periods. 

seen in the right figure of Fig. 1. In order to determine how well the averaged 
equation (36) captures the differences between the purely diffusive case and the 
case with both time-dependent advection and diffusion, the difference of the 
purely viscous case and the solution of (1), hence c—c^^^^^ and (b) the difference 
of the purely viscous case and the approximation c^^^^ constructed using (36), 
hence plotted in Fig. 2. The first order approximation accurately 

captures the overall dynamics of the full, time-dependent problem. 




20 40 60 80 TOO 120 20 40 60 BO 100 120 



Fig. 2. Comparison of the prediction of the averaged equation (36) to the evolution of 
the time-dependent equation (1). We plot the difference between the purely viscous 
solution and the solution incorporating the effects of the time-dependent advection 
field. The figure on the left shows c — where c is the solution of (1). The figure on 
the right shows c^^^^ — c*^^'^-* , where c^^^^ has been found using the averaged equation 
(36). The averaged equation is obviously able to capture the leading order impact 
of the time-dependent velocity field. 

To quantify the accuracy of the approximation, we consider the time-evolution 
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of the canonical L^-norm of the differences by defining 



c — c 



{av)\ 



/lR2 


c — 


C 


(av) 


^ dx dy 


/lR2 


C 


^ dx dy 



1/2 



(43) 



Figure 3a shows that this error is small and approximately constant over 
the first 10 periods whereas the corresponding error between c and c^^'*^-* is 
comparatively large and growing exponentially in time. Figure 3b indicates 
that, at least for the case where f(t) = sin(t), the solutions of the averaged 
equation converge to those of the full equation faster than e. 




e = T x/ L-" 



Fig. 3. (a) Difference norms for the case f{t) = sin(t) and e = 0.005. (b) Average 
error between the first order averaged equation and the full solution versus e for 
f{t) = sin(t). The error scales like e^'^. 

We examine the role of f{t) in the averaged dynamic by setting f{t) = cos(t). 
This choice implies that (F) = 0, leading to near degeneracy of the first order 
corrections. In this case, the time independent advection terms produced by 
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the averaging procedure do not contribute to 'twisting' the scalar evolution. 
The symmetry breaking terms in the averaged diffusivity tensor are identi- 
cally zero at first order. The effect of this near-degeneracy for (F) = is 
clearly seen in the comparison of panels (a) and (b) in Fig. 4. A comparison 
of the difference norms, shown in Fig. 5(a), indicates that the first order aver- 
aged equation is only a marginal improvement on the purely viscous solution 
at short times. However, as shown in Fig. 5(b), solutions to the first order 
averaged equations continue to converge to the true solution with decreas- 
ing e although the convergence rate, ~ e'-\ is considerably slower than that 
observed when f(t) = sin(t). 



(a) 





(d) 



Fig. 4. Solutions for the case f{t) = cos(t) and e = 0.010 after 10 periods, (a) Full 
solution, (b) first order averaged solution, (c) second order averaged solution, (d) 
difference between first and second order averaged. 

For (F) = 0, second order contributions are clearly important. Referring back 
to (40), this situation also leads to a relatively simple form for the second 
order expression. The coefficient in front of £i vanishes, and in this particular 
case for f(t) = cos(t), (40) simplifies to 



Vt = e{AV+-C2V]+e' 



[A,C,]V-^-[C,,C2]V^ 



(44) 



13 



Evaluating this equation explicitly we find 



Vt = e{AV+{uj')Xe)+^' 




) 



(45) 



As shown in the Appendix, this leads to 0{e) symmetry breaking contribu- 
tions to both the average advection terms and diffusivity tensor as well as 
contributions in the form of higher, third order, spatial derivatives. Numeri- 
cally, such terms are easily computed spectrally. An example of a second order 
solution is shown in panel (c) of Fig. 4. Obviously, the second order solution 
is a clear improvement on the first order approximation and the difference be- 
tween the two, shown in panel (d), demonstrates the restoration of advective 
twist at higher order. Figure 5 indicates both the increase in accuracy and the 
expected increase in convergence rate for the second order approximation. 



7 Discussion 

We have proposed a scheme for formally transforming the advection diffusion 
equation, in the limit of small diffusivity, into a form suitable for averaging. 
We have given an explicit means of averaging the transformed equation and 
proven the convergence of solutions of the averaged, time-independent ap- 
proximation to the full dynamics. Throughout, however, we deal only with 
a restricted class of advecting fields, namely those which are zero-mean and 
possess time-periodic stream lines. While such fields are inherently integrable, 
and hence explicitly non-chaotic, the results presented are of interest in the 
study of 'strange eigenmodes' of the advection diffusion equation in the Batch- 
elor regime. First, the emergence of strange eigenmodes is independent of the 
integrability or non-integrability of the underlying flow [5], and indeed, the 
simple example considered here produces non-trivial periodic patterns. Sec- 
ondly, the results shown for even these extremely simple cases point to the 
delicate relationship between the non-linearity of the conservative advection 
operator and small diffusivity. The analysis points to fundamental differences 
between the dynamics of continuous time flows and discrete time maps. In- 
deed, for the periodically modulated vortex considered, the averaged dynamics 
depends strongly on the phase of the single frequency periodic modulation, a 
fact completely lost when considering the Poincare map of the flow which is 
simply the identity. 

In the case of a mean-free advection field with periodic stream lines, the trans- 
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e = T k/ 

Fig. 5. (a) Difference norms for tlie case f{t) = cos(t) and e = 0.005 for both first 
and second order averaging, (b) Average error between the averaged equation and 
the full solution versus e for f{t) = cos(i). The first order solution is shown by the 
solid curve, the second order by the dashed curve. For comparison, the first order 
results for f{t) = sin(i) are shown in the light dot-dashed line. 



form to action-angle variables of the Lagrangian flow was found to be the ap- 
propriate transformation for deriving accurate time-averaged dynamics in the 
limit of small diffusivity. Finding transformations with similar properties for 
other classes of advection fields will likely provide a means for understanding 
the interplay between advection and small diffusion. 
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A The averaged equation in Cartesian coordinates 

In order to convert (36) back to Cartesian coordinates, we note first that it 
follows directly from the transformation rules to polar coordinates that 



which, after straightforward calculations, yield the following transformation 
rules for the higher operators occurring in (36): 

dede = -xd^ - ydy 

+y'^da:dx - 2xydxdy + x^dydy 
rdedr = xdy - yd^ 

-xydA + {x"^ - y^)dxdy + xydydy 

Rescaling time as r = et, we find for (36) then 



rdr = xdx + yd^ 



dg = xdy - yd, 



'X 



Vr + (E- V)V = K : VVV 



(Al) 



where the advection vector S and is given as 



'^2{x,y) 



^i{x,y) 




and the diffusion matrix K becomes 



Kn = -2-xy{F) + {F'){u;')Y 
r 

K,^ = K^i = -{F){x'-y') - {F')iu;rxy 
r 

K22 = 2-xy{F) + {F'){uj'fx' 



16 



This together with the formulas 



ijj = 9, \- UJ 



(a2 + r2)2' r {r^ + a^^ 

gives directly the coefficients of (36) in Cartesian coordinates. Note that all 
operators give a contribution to the averaged advection field in Cartesian 
coordinates. 

For the second-order, the corresponding terms are more complicated. The 
explicit transformations of the operators are given by 



r'ded^, = 2xydl + 2{x'' - y')d,dy - 2xydl 

+xy^dl + (x^^ - 2xy^)dldy + {2x^y - y^)d,dl - x^ydl 
dl = -xdy + yd^ 

+3xydl + 3{y^ - x'^)dxdy - 3xydy 
Wd^ - Sx^yd.d^ + 3xy^dldy - y^d^ 

Note in particular that the ffist term in the transformation of dg corresponds 
to -de, hence is one example of a term that will have a "twisting" effect on 
the solution and that can lead to symmetry breaking. We can now rewrite 
(45) in Cartesian coordinates (rescaling time r = et as before) and obtain: 

K + (S ■ V) F = KuV,, + 2KuV,y + K22Vyy 

+ MuiV^XX + MluV^^y + Ml22V^yy + M222Vyyy (A.2) 

where the coefficients of the advection field S are given by 

^i = hx + e[{h~h)y + f2^ 
S2 = /5Z/ + e((/3-/i)x-/2^), 
and the diffusion matrix K results in 

K,, = l + hy' + e (Shxy - /^^ - 2/4^) 
K22 = l + hx' + e (-3hxy + /^^ + 2/4^) 

2K,2 = -2hxy + e (/2^^ + 3/3(1/' - x') + 2/4^^) . 

As for the tensor M yielding higher-order contributions, we find 
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2 

and the coefficients are explicitly written as 



, , // , /" ^4 /I ^2^2 _|_ „4 

, , /a;" 8 9r^ - 14a + 



y r y r (r^ + a^)^ 

/3 = 4- - 2{u'fu" = - . + 16 

- 2a^ 



/4 = 4u;" = 4 



(r2 + a?f 



1 



2' ' (a2 + r2)4 
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